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Abstract 

°Li atoms in the three lowest hyperfine states display universal properties when the S-wave 
scattering length between each pair of states is large. Recent experiments reported four pronounced 
features arising from Efimov physics in the atom-dimer relaxation rate, namely two resonances and 
two local minima. We use the universal effective field theory to calculate the atom-dimer relaxation 
rate at zero temperature. Our results describe the four features qualitatively and imply there is a 
hidden local minimum. In the vicinity of the resonance at 685 G, we perform a finite temperature 
^ calculation which improves the agreement of theory and experiment. We conclude that finite 

temperature effects cannot be neglected in the analysis of the experimental data. 

PACS numbers: 31.15.-p,34.50.-s, 67.85.Lm, 03.75.Ss 
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I. INTRODUCTION 



Few-body systems with large two-body scattering length a are of great interest because 
they have universal properties insensitive to the details of the interaction at short distances. 
The simplest example is a universal two-body bound state with binding energy h 2 /(ma 2 ) 
when a > 0. The corrections to this formula are suppressed by £/a where t is the natural 
low-energy length scale. For alkali atoms, i is given by the van der Waals length £ v dw which 
quantifies the range of the interaction. Examples of universal two-body states can be found 
in different branches of physics such as atomic physics, nuclear physics, and high-energy 
physics [U |2]. Universal properties also exist for more than two particles. Systems of three 
identical bosons have an infinite number of geometrically spaced low-energy bound states 
(Efimov trimers) with an accumulation point at zero energy when the scattering length is 
taken to infinity and the range of the interaction is taken to zero [3] . In this so-called unitary 
limit, the binding energies of two successive trimers are related by a multiplicative factor of 
approximately 22. 7 2 ~ 515. This remarkable feature is the consequence of a discrete scaling 
symmetry with a scaling factor e n ^ s ° ~ 22.7 where so ~ 1.00624 for identical bosons. For 
more complicated systems, the scaling factor depends on the mass ratios and the spin of 
the components. We will refer to phenomena that result from the discrete scale invariance 
in the three-body sector as Efimov physics [J]. Corrections to Efimov physics in a process 
with typical wave number k are suppressed by tja and by ka. 

Efimov trimers in ultracold atomic gases can be observed through enhanced loss rates in 
three-body recombination and atom-dimer relaxation processes [5H7J . The first evidence for 
an Efimov trimer in a system of bosonic 133 Cs atoms was presented by Kraemer et al. [8]. In 
a subsequent experiment with a mixture of 133 Cs atoms and dimers, Knoop et al. observed 
a resonance in the loss of atoms and dimers [9] which can be explained by an Efimov trimer 
crossing the atom-dimer threshold [10J. Efimov trimers have also been observed using other 
types of bosonic atoms. Zaccanti et al. [11] found evidence for two adjacent Efimov trimers 
in a gas of ultracold 39 K atoms and Gross et al. [12] measured two Efimov features across 
a Feshbach resonance and confirmed the universal relation between the features with 7 Li 
atoms. Attached to each Efimov trimer are two universal tetramers which can be observed 
in four-body loss processes [T3HT5] . Ferlaino et al. [H] observed signatures of two such 
tetramers close to an Efimov trimer with 133 Cs atoms. In 7 Li atoms, Pollack et al. [17] have 
seen two sets of two tetramers that are close to the corresponding Efimov trimers. 

Efimov physics can also occur in fermionic systems with three or more spin states [3j. 
Several features associated with Efimov trimers in 6 Li atoms have been reported recently. 
Two three-body recombination features were observed by Ottenstein et al. [18] and Huckans 
et al. [19]. These features were analyzed theoretically in Refs. [201422] and traced back 
to an Efimov trimer close to the three-atom threshold [2TH2U1J . Another recombination 
feature was discovered by Williams et al. [24J and identified as a second trimer close to the 
threshold. Braaten et al. [22] analyzed the energy spectrum of Efimov trimers and predicted 
the crossings of the trimers with the atom-dimer threshold. These crossings imply two 
resonances in atom-dimer scattering. The resonances were identified in recent experiments 
by Lompe et al. [25] and Nakajima et al. [26J. Furthermore, two local minima in the atom- 
dimer relaxation rate were reported in Ref. [25] . 

In this paper, we study Efimov physics in atom-dimer scattering of 6 Li atoms in the three 
lowest hyperfine states and perform a complete zero-range calculation at zero temperature. 
The size of finite range corrections to our results is estimated. We identify the two resonances 
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and the two minima in the atom-dimer relaxation rate. Quantitatively, our zero temperature 
results show deviations from the data. An approximate finite temperature calculation near 
one of the resonances improves the description of the data considerably and shows that finite 
temperature effects account for about 25% of the discrepancy. The paper is organized as 
follows: In Sec. Inland Sec. Ill we explain our theoretical framework. The numerical results 



at zero temperature are displayed and compared with the experimental data in Sec. IV In 



Sec. [Vj we carry out an approximate finite temperature calculation near the resonance at 
B = 685 G. We summarize our results and conclude in Sec. IVD 



II. DIMER RELAXATION 

In this section we discuss the losses of atoms and dimers through inelastic scattering 
processes and present expressions for the relaxation rate constants. 

Before proceeding with our discussion of atom-dimer scattering, we explain our notation 
and terminology. We label an atom in one of the three hyperfine states of the 6 Li atoms 
with an index i where i — 1, 2 or 3. The S-wave scattering length between atoms in states i 
and j is denoted either as or as where k ^ i ^ j. Two atoms in the same state can not 
scatter in an S-wave because of the Pauli principle. If the scattering length is positive 
and much larger than the van der Waals length £ v dw ~ 65ao, the atoms i and j can form a 
dimer with binding energy h 2 /(ma 2 ^). We call this dimer the (shallow) r/'-dimer or simply 
a shallow dimer. Shallow dimers have to be distinguished from deep dimers with binding 
energy of order h 2 j (m^dw) or l ar g er - 

In a gas of atoms i and jfc-dimers, the atoms and dimers can undergo inelastic collisions 
into atoms and deeply bound dimers with a binding energy larger than that of the jfc-dimer. 
This inelastic process is called dimer relaxation. The difference in the binding energies of 
the initial and final state dimers is released as kinetic energy and the atom and dimer in 
the final state recoil from each other. If their kinetic energies are larger than the trapping 
potential, they escape the trap. The loss rate for the number density rii of atoms i and 
number density rijk of jfc-dimers is 

dl ni = dl njk = ~^ irLin i k ' w 

where the coefficient is the relaxation rate constant for the j/c-dimer and atom i. 

In the case of identical bosons, dimer relaxation is possible only if the final state consists 
of an atom and a deep dimer. However, in the three-fermion system relaxation can also 
occur into into shallow dimers. For example, in a scattering process of an atom i and a jk- 
dimer, relaxation can proceed into the r/'-dimer provided the binding energy of the r?'-dimer 
is larger than that of the jfc-dimer. Therefore, the total rate is the sum of all relaxation 
rates into atoms plus shallow dimers and atoms plus deep dimers 

& = E^V^ deep > ( 2 ) 

where the index i implies atom i plus jfc-dimer in the initial state and the j implies atom 
j plus ifc-dimer in the final state. For brevity, we refer to as the rate into the shallow 

dimer or the rate into the (final) ifc-dimer and to (3f eep as the rate into the deep dimers. 



3 



The relaxation rate can be calculated from the T-matrix element for atom-dimer scat- 
tering 7if D (k,p; E) where k and p are the relative wave numbers of the atom and dimer in 
initial and final state, respectively, and E is the total energy. By using the optical theorem, 
we can calculate the total rate for atom-dimer scattering which is the sum of elastic and 
inelastic rates. However, in the low energy limit k — > 0, the elastic rate scales as k because 
7^f D (A;, k; E) is constant and the two-body phase space gives one power of k [I]. The elastic 
rate therefore vanishes at zero energy and the optical theorem gives the total relaxation rate 

Of, 

A = ^7-^(0,0,-1/(77^)). (3) 

The rate into the shallow dimer (3f\j is determined by the square of the T-matrix element 
multiplied by the two-body phase space: 



Ptj = |^ \r4°(0,p, -1/W))| 2 Oifii - aj), ( I) 



where the wave number p = (2/\/3)y a j 2 ~ a i 2 an d the 8- function is inserted because the 
relaxation is allowed only when the initial dimer binding energy is smaller than the final 
one. The prefactor 2p/(37r) comes from the two-body phase space integral. From Eq. (|4]), 
one can deduce that /3f^ vanishes as ^a-i — aj near the crossing of the ik- and j/c-dimers 
where a» approaches aj. 

The relaxation rate into deep dimers (3f eep could also be calculated by using the T-matrix 
in a similar way to Eq. Q if the theory described deep dimers explicitly. Alternatively, the 
effects of deep dimers can be taken into account indirectly by using an analytic continuation 
of the three-body parameter into the complex plane as introduced in Ref. [27J. Then, the 
partial rate fS^ eep can be obtained by using the relation in Eq. 

D'Incao and Esry have calculated the scattering length dependence of ultracold three- 
body collisions near overlapping Feshbach resonances for a variety of cases [28]. For the 
relaxation rate into the shallow dimer /S£_^-, they find that it scales like a 2 jai times a log- 
periodic function of aj if ^> aj while at is non-resonant. The rate into deep dimers 
/3f eep scales like a 2 /ai times a prefactor that is a constant for positive aj and a log-periodic 
function of aj for negative aj. 



III. STM EQUATION 

The Skorniakov-Ter-Martirosian (STM) equation [29] is an integral equation that de- 
scribes three-atom scattering interacting through zero-range interactions. In this section, 
we discuss the STM equation and relate the T-matrix element to the amplitude Aij 
that is the solution to the STM equation. 

We consider only the S-wave contribution and assume that higher partial wave contribu- 
tions are suppressed at low temperature. For the three-fermion system, the STM equation 
forms 9 coupled equations for the amplitudes Aij(k,p; E) [201 E3]. For non-zero energy, the 
equation is given by 

Aij(k,p;E) = (l-5 tJ )Q(k,p;E) 

+ Ij2(l-S kj ) [ dqq 2 Q(q,p;E)D k (q;E)A lk (k,q;E), (5) 
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where Aij(k,p; E) is the amplitude for an atom % and a complementary pair of atoms to 
scatter into an atom j and a complementary pair and A is an ultraviolet cutoff. The function 
Q(k,p; E) and the 2-atom propagator Dk(q; E) are given by 



Q{k,p; E) 
D k (q; E) 



2kp 



In 



k 2 + kp + p 2 — mE — ie 
k 2 — kp + p 2 — mE — ie 
1 



-l/a k + \l\q 2 - mE - ie 



(6) 
(7) 



The solutions of the STM equation ^ depend log-periodically on A with a discrete 
scaling factor e n ^ s ° ~ 22.7, where Sq ~ 1.00624. The dependence on the arbitrary cutoff A 
can be eliminated in favor of a physical 3-body parameter such as the binding wave number 
of an Effmov trimer in the unitary limit. For convenience, we choose to work directly with 
the wave number cutoff A in our calculations. If deep dimers are present, the trimer has 
a finite width that allows it to decay into an atom and a deep dimer. The effects of deep 
dimers can be taken into account by analytically continuing the cutoff A into the complex 
plane [27] 

A_>A e Ws 0) ( 8 ) 

where 77* is a width parameter associated with the effects of the deep dimers. A and 77* 
are determined from experimental measurements of three-body recombination in the 6 Li 
system [[24] and their numerical values are A = 456 ag 1 and 77* = 0.016 [23] where ao is the 
Bohr radius. 

In order to obtain the T-matrix element, the amplitude Aij must be multiplied with the 
dimer- wavefunction renormalization factors -J ZiZj^ where Z\ = 87r/aj, 



(9) 



By solving the STM equation ([5| and using the relation in Eq. ([9]), we can calculate the 
relaxation rates in Eqs. ^ and ^4]). 



IV. ZERO TEMPERATURE RESULTS 

In this section we present our numerical results for the dimer relaxation rate constants 
at zero temperature and compare them with recent measurements in Refs. [251 126] . We 
solve the STM equation in Eq. ^ numerically with 5 input parameters: the three pair 
scattering lengths a^, % = 1,2,3, the cutoff A and the width parameter 77*. We use the 
values A = 456 1 and 77* = 0.016 that have been determined from measurements of the 
three-body recombination rate in Ref. [24] and have been used in Ref. [23] to predict the 
atom-dimer relaxation rate. 

Fig. [T] shows the scattering lengths for the 3 lowest hyperfine state of 6 Li as a function of 
the magnetic field. The atom-dimer relaxation rate is non-zero in the region of positive 
scattering length cij. The upper limit of the region is set by the Feshbach resonances at 
811 G, 690 G, and 834 G for a\, (12, and a 3 , respectively. If the scattering length is much 
larger than the van der Waals length £ v dw ~ 65a , the universal theory is valid. We denote 
this region as the universal region. Corrections due to the finite range of the interaction 
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FIG. 1: (Color online) The scattering lengths for the three lowest hyperfine states of 6 Li in units 
of 10 3 ao as a function of the magnetic field B from [30]: an = a 2 (long-dashed), avi = «3 (short- 
dashed), and = a% (dashed-dotted) . The vertical lines at 690 G, 811 G, and 834 G indicate the 
Feshbach resonances for 013, 023, and 012, respectively. 



should be small in the universal region. In practice, we apply the universal zero-range theory 
when all scattering lengths are at least two times larger than £ v dw corresponding to magnetic 
fields B > 608 G. The expected error due to finite range corrections is £ v dw/a an d therefore 
smaller than 50% in this region. 

In Ref. [23J, two crossings of Efimov trimers with the atom-dimer threshold have been 
predicted. Both are located at the l(23)-threshold: at B* = 672 G and at B'^ « 597 G. Here 
the index 1 denotes the atom and the index (23) the dimer. Since the atom-dimer relaxation 
is resonant when the trimer appears near the threshold, two resonances are expected in the 
rate f3\ near B* and B'^. The resonance at B* is well in the universal region where all 
scattering lengths are much larger than £ v aw while B^ is slightly outside. Therefore, the 
resonance position should be accurately determined with corrections of order £ v dw/ a 23 ~ 
10% where the value of 023 at the resonance was used. The position is outside the 
universal region and can receive large non-universal corrections of order 100%. These error 
estimates are accurate up to a prefactor of order one. The exact value of this prefactor can 
only be obtained from an explicit calculation of the range corrections. Note also that these 
percentage errors apply to the positions in terms of the scattering length. To obtain the 
errors for the corresponding magnetic field they have to be converted using Fig. [TJ 

In Fig. [2j we show our numerical results for /3\ and compare them with the recent measure- 
ments of Refs. [25j [26] . We give the full relaxation rate as well the individual contributions 
from shallow and deep dimers. In the magnetic field region from 590 G to 730 G, only 
the relaxation into deep dimers contributes to /3\ since the energy of the 23-dimer is larger 
than the energy of the other shallow dimers. At a magnetic field of 730 G, the 23-dimer 
crosses the 12-dimer and the relaxation channel into 12-dimers opens up. As the crossing 
is approached from above, the rate vanishes as \fa-\_ — 03 in agreement with the analytical 
result from Sec. [XTJ Between 790 and 810 G, the condition a\ 3> with a 2 non-resonant is 
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FIG. 2: (Color online) The relaxation rate constant for the 23-dimer and atom 1 as a function of the 
magnetic field B. The squares and circles are data points from Ref. [25] and Ref. |26j . respectively. 
The curves are our results for the total rate f3\ (solid line), the partial rate into atom 3 and the 
12-dimer /3f t \ K3 (dashed line), and the rate into an atom and a deep dimer /3^ eep (dashed-dotted 
line). The light solid line gives the total rate /3i for the parameters obtained in Ref. [25] while the 
vertical line marks the boundary of the universal region. 

satisfied approximately. The rate into the shallow dimer in this region scales approximately 
as a|/ai in agreement with the prediction of D'Incao and Esry [28]. At 811 G, the 23-dimer 
disappears through the three-atom threshold and the relaxation rate vanishes. 

Our results show two resonances at £>* = 672 G and B^ 597 G. There is also a dramatic 
change in the relaxation rate at 730 G because the relaxation channel into 12-dimers opens 
and the corresponding rate into the 12-dimer increases rapidly. Our results describe the 
resonances in the experimental data qualitatively. The second resonance has been measured 
at -B* xp = 685 G in Ref. [251 ES]- This value is 13 G away from the theoretical prediction 
£>* = 672 G. In terms of the scattering length, this corresponds to c^-Bf^Vc^-B*) = 
1076/835 ~ 1.3, leading to a 30% shift in the resonance position. This shift is a factor three 
larger than the naive error estimate of 10%. Taking into account the unknown prefactor of 
order one in the estimate, however, the two values are consistent. Except near the resonance, 
the experimental data are generally above our results. 

Lompe et al. [25] analyzed their data near the resonance at using an approximate 
analytic expression from [23] and extracted the resonance position and the width parameter 
r/*. They found the value r/* = 0.34 which is more than an order of magnitude larger than the 
value 0.016 extracted from three-body recombination. Moreover, the normalization of the 
relaxation rate was adjusted to describe the data. The light solid curve in Fig. [2] gives our 
universal result for the parameters A = 329 ag 1 (which reproduces the resonance position 
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FIG. 3: (Color online) The relaxation rate constant for the 13-dimer and atom 2 as a function of 
the magnetic field B. The squares are data points from Ref. [25] ■ The curves are our results for the 
total rate fa (solid line) and the partial rates into an atom 3 and a 12-dimer /3|^. 3 (dashed line), 
into an atom 1 and a 23-dimer (3^ l (dashed-dotted line), and into an atom and a deep dimer fl^^ 
(dotted line). The light solid line gives the total rate /3i for the parameters obtained in Ref. |25j 
while the vertical line marks the boundary of the universal region. 

^exp _ ggg q-j an( ^ ^ _ q 34 These parameters give a much better description of the 
data in Refs. [25j EH] but are generally a factor 2-3 above the experimental data. However, 
one should keep in mind that the experimental data are only a factor three lower than the 
unitarity bound and finite temperature effects are likely important. We will come back to 
this issue in the next section. 

The reason for the considerably larger value of 7/* extracted in [23] compared to the value 
from recombination data is not understood. However, we note that a similar discrepancy 
between the values of 77* from atom-dimer relaxation and three-body recombination occurs 
in the bosonic system of 133 Cs atoms [HI ITU] . 

Nakajima et al. [26J performed a numerical analysis for /3i based on the universal theory 
and on a two-channel model. Their results obtained with the universal theory agree with 
our calculation. Within the two-channel model they derived energy-dependent scattering 
lengths which introduce non-universal effects in the two-body amplitudes. Because this 
model cannot resolve the discrepancy between the universal results and the measurements 
they concluded that the three-body parameters A and r/* depend on the magnetic field. 

In Fig. [3j we show the experimental data from Ref. [23] and our numerical results for /3 2 . 
The total rates, both of the data and the numerical results, show no pronounced structure. 
However, there is a local minimum in the partial rate near 600 G that is outside 

the universal region. As discussed in Ref. [28], this minimum is the effect of destructive 
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FIG. 4: (Color online) The relaxation rate constant for the 12-dimer and atom 3 as a function of 
the magnetic field B. The squares are data points from Ref. [25]. The curves are our results for 
the total rate /?3 (solid line), the partial rate into atom 1 and a 23-dimer (dashed line), and 

the rate into an atom and a deep dimer /3g 6ep (dashed-dotted line). The light solid line gives the 
total rate f3\ for the parameters obtained in Ref. [25] while the vertical line marks the boundary 
of the universal region. 

interference between different recombination channels. This interference pattern does not 
exist in systems with identical bosons where relaxation can occur only into deep dimers. 
In the total rate, the interference is hidden by the dominant process /3|^ >3 . The light solid 
curve again gives our universal result for the parameters A = 329 Oq 1 and 77* = 0.34 obtained 
in Ref. [25J. The difference between the two parameter sets is very small for /3 2 , but the 
alternative set gives a slightly better description of the data. Between 670 and 690 G, a 2 is 
much larger than a\ and 03 which are approximately equal. The scaling of the calculated 
rates into the shallow dimers is consistent with the prediction of D'Incao and Esry [28J: 
$£> 3 ~ a l/ a 2 and /?^ 3 ~ a\/a 2 . 

Fig. [4] shows our results and the recent measurement [25] for (3%. Between 590 G and 
730 G, the energy of the 12-dimer is smaller than the energy of the 23-dimer and the 
relaxation channel into 23-dimers is open. After the 23-dimer crosses the 12-dimer at 730 G 
this relaxation channel is closed and only relaxation into deep dimers is possible. As the 
crossing is approached, the rate vanishes as i/a 3 — aq in agreement with the analytical 
result from Sec. [n] Between 832 and 834 G, a x is negative and the condition a 3 ^> |ai| 
with a 2 non-resonant is satisified approximately. The rate into deep dimers scales in this 
region approximately as af/a 3 consistent with the prediction by D'Incao and Esry. Two 
interference minima have been observed at 610 G and 695 G [25J while our results show two 
minima at 600 G and 715 G. Above 700 G, the data are larger than our results by more 
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than a factor 10. Using the alternative parameters A = 329 1 and 77* = 0.34 obtained in 
Ref . [25J , we again find a better agreement with the data. With these parameters the second 
minimum in the rate into the shallow dimer disappears beyond 730 G. Hence, the partial 
rate decreases monotonically and vanishes near 730 G. For smaller 77*, the position of the 
minimum in the total rate is around 730 G. As 77* increases, the position remains almost the 
same and the depth of the minimum becomes shallower. The minimum is not visible when 
77* = 0.34 because the rate into deep dimers is much larger than the rate into the shallow 
dimer. Therefore, with these parameters that were fit to data for /3i the position of the 
second minimum in the data for /3 3 cannot be explained correctly. 

If the rate into 23-dimers could be separated experimentally from the total rate, it would 
clearly determine the positions of the local minima. This could be achieved by tuning the 
depth of trapping potential such that it is much larger than the kinetic energies of atoms and 
23-dimers in the final state but much smaller than the energies of atoms and deep dimers in 
the final state. The kinetic energies of an atom and a deep dimer in the final state could be 
estimated from the binding energy of the deep dimers. Their energies would be of the order 
of the van der Waals energy or larger: E v $w/h ~ 154 MHz, where -E v dw — ^VO^vdw)- 1 
The kinetic energies of an atom and a 23-dimer are given by the difference in binding energies 
between the 12-dimer and the 23-dimer: %/ {^l-nm)^^ — a^ 2 2 ) is about 1 MHz at 650 G and 
vanishes at 730 G. This way one may be able to measure the rate into deep dimers separately 
and to extract the rate into the 23-dimer. 



V. FINITE TEMPERATURE RESULTS 

The results from the previous section suggest that finite temperature effects may play an 
important role in understanding the atom-dimer relaxation data from Refs. j25j [26] . A full 
finite temperature calculation of the 6 Li system is beyond the scope of this work. Therefore, 
we perform an approximate calculation of the relaxation rate Pi near the resonance at 
BfP = 685 G. 

We start from the approximate expression for the scattering length between an atom 1 
and a 23-dimer near the resonance at that was extracted from a calculation of the trimer 
binding energy in |23j : 

01(23) ~ (Cicot[s ln^/a*) + + C 2 ) a 23 , (10) 

where a* = 023(5*) and the coefficients are C\ = 0.67 and C 2 = 0.65. Using the scattering 
length approximation for the S-wave atom-dimer scattering amplitude in the l(23)-channel, 

/i(23)(fc) = [-l/oi(23) (11) 

we can calculate the inelastic scattering cross section. At low temperatures, the contributions 
from higher partial waves can be neglected. We subtract the elastic cross section from the 
total cross section obtained via the optical theorem as in [31J and find for the inelastic cross 
section: 

(inelastic) /, s _ 4?T ~ ima l(23) ^ „x 

CTl(23) 1 ] ~ k 1 - 2Adma 1(23) + fc 2 |a 1(23) r 1 ] 



A convenient conversion constant for 6 Li atoms is given by % /(mag) = 6OO/1 GHz= 28.8fce K. 
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FIG. 5: (Color online) The relaxation rate constant /3\ for a 23-dimer and an atom 1 as a function 
of a/a* near the resonance at 5* xp = 685 G (a* = 1076 ao). Squares and circles are data points 
from Ref. [53] and Ref. [2B], respectively. The solid (dashed) curves correspond to 77* = 0.016 
(?7* = 0.34). The upper curves give the zero temperature result, while the lower curves give the 
finite temperature result for T = 60 ± 15 nK with the shaded area indicating the temperature 
uncertainty. 

The total dimer-relaxation rate 0i can then be calculated by taking a Boltzmann thermal 
average of the inelastic reaction rate v re i %1 (fc) where v re \ = 3hk/(2m) is the relative 
velocity of the atom and dimer in the initial state. This leads to the expression 

ftCO = ^(ka^\k)) = I^Vf [ k^dk ka^\k) e~&, (13) 

where At = \J 2ith 2 / (rnksT) is the thermal de Broglie wavelength of the atoms. 

Our results for the atom-dimer relaxation rate constant near the resonance at -B* xp = 
685 G are shown in Fig. [5j The results for the width parameters 77* = 0.016 and 77* = 0.34 
are given by the solid and dashed curves, respectively. In each case the upper curves give 
the zero temperature result while the lower curves give the finite temperature result for 
T = 60 ± 15 nK. Here, the shaded area indicates the uncertainty from the temperature. The 
value of a* = 1076 a has been fixed to reproduce the resonance position at Bf p = 685 G. 
For 77* = 0.016, the finite temperature effects decrease the height of the peak by an order 
of magnitude but the predicted resonance shape is still much narrower than the data. The 
finite temperature effects are much less severe for 77* = 0.34. They only lead to a reduction 
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of Pi by about 25 % but clearly improve the description of the data. We conclude that 
finite temperature effects can not resolve the question of the different values for 77* in the 
three-body recombination and dimer-relaxation data. Moreover, while finite temperature 
effects are important, a qualitative description of the data at T 60 nK can be already 
achieved with a zero temperature calculation. 

VI. SUMMARY AND OUTLOOK 

In this work, we have studied Efimov physics in atom-dimer relaxation of 6 Li atoms in 
the three lowest hyperfine states using the universal zero-range theory. Two resonances were 
observed at magnetic fields 603 G and 685 G in the relaxation rate f3\ in recent experiments 
[251 [26] . These resonances are consequences of two Efimov trimers close to the atom-dimer 
threshold. Their positions have been predicted by Braaten et al. [23]. The measured position 
of the resonance at 685 G [25] [26] , which is well within the universal region, is larger than 
the prediction by about 30%. This is consistent with an error of order 10% error due to 
effective range corrections. However, the value 77* = 0.016 extracted from the three-body 
recombination data is not able to describe the atom-dimer relaxation data which require 
the larger value 77* = 0.34 [22]. The reason for the larger value of 77* in dimer relaxation 
compared to the value from recombination data is not understood. However, we note that 
a similar discrepancy between the values of 77* from atom-dimer relaxation and three-body 
recombination occurs in the bosonic systems of 133 Cs atoms pi HH] ■ 

Using the value 77* = 0.34, our zero temperature calculation is able to describe the data 
qualitatively. In the vicinity of this resonance, we have also performed an approximate finite 
temperature calculation and find sizable temperature effects that can suppress the relaxation 
rate by an order of magnitude if 77* = 0.016. For the larger value 77* = 0.34, however, these 
effects lead to a moderate suppression of about 25%, such that zero temperature results are 
useful as a first approximation. 

In Ref. [25], also two local minima at 610 G and 695 G are discovered in the rate ^3. Those 
minima can be associated with destructive interference between different recombination 
channels [28] . Our numerical results show that the partial rate into 23-dimers is responsible 
for the minima but the positions of the minima are displaced by -10 G and +20 G from the 
measurements. These displacements correspond to a 30% shift in terms of the scattering 
length compared to our predictions. Since the scattering length, a 2 3, is about a factor two 
larger in this magnetic field region than around the resonance in j3i, one would expect smaller 
corrections here. The observed shifts indicate that the destructive interference leading to 
the minima might be dominated by wave numbers k larger than 1/023 such that corrections 
of order k£ v dw are important. Moreover, finite temperature effects could fill the minima in 
an asymmetric fashion. The total rate /3 2 shows no structure. However, the partial rate into 
the 23-dimer shows a local minimum near 600 G. In the total rate this feature is hidden by 
the dominant rate into the 12-dimer. 

In order to better understand the discrepancy between the values of 77* and the reso- 
nance positions extracted from three-body recombination and atom-dimer relaxation two 
important improvements are required in a future analysis. First, a full finite temperature 
calculation of atom-dimer relaxation should be carried out. This requires calculating atom- 
dimer scattering above the dimer breakup thresholds and in higher partial waves. Due to 
the different pair scattering lengths in the three channels, such a calculation is considerably 
more complex than in the case of identical bosons [32J. However, it will allow to better 
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distinguish effects from the resonance width parameter 77* and from the finite temperature 
which can be partially traded for each other [Till 132] . Second, an analysis of the effective 
range corrections should be performed in order to describe the observed shifts in the reso- 
nance positions quantitatively. A similar analysis for systems of identical bosons was carried 
out in Refs. [331 EU- 
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